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ABSTRACT 



We used Early Release Science (ERS) observations taken with the Wide Field Camera 3 (WFC3) in the GOODS-S field to study the 
galaxy stellar mass function (GSMF) at 0.6 < z < 4.5. Deep WFC3 near-IR data (for Y as faint as 27.3, J and H as faint as 27.4 
AB mag at 5a), as well as deep Ks (as faint as 25.5 at 5a) Hawk-I band data, provide an exquisite data set with which determine in 
an unprecedented way the low-mass end of the GSMF, allowing an accurate probe of masses as low as M, = 7.6 • 10 9 M o at z ~ 3. 
Although the area used is relatively small (~ 33 arcmin 2 ), we found generally good agreement with previous studies on the entire mass 
range. Our results show that the slope of the faint-end increases with redshift, from a = - 1.44 ± 0.03 at z ~ 0.8 to a = -1.86 ± 0.16 at 
z ~ 3, although indications exist that it does not steepen further between z ~ 3 and z ~ 4. This result is insensitive to any uncertainty 
in the M* parameter. The steepness of the GSMF faint-end solves the well-known disagreement between the stellar mass density 
(SMD) and the integrated star formation history at z > 2. However, we confirm the that there appears to be an excess of integrated star 
formation with respect to the SMD at z < 2, by a factor of ~ 2 — 3. Our comparison of the observations with theoretical predictions 
shows that the models forecast a greater abundance of low mass galaxies, at least up to z ~ 3, as well as a dearth of massive galaxies 
at z ~ 4 with respect to the data, and that the predicted SMD is generally overestimated at z < 2. 
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1. Introduction 

Understanding the assembly of stellar mass in galaxies is a fun- 
damental step towards a description of galaxy evolution. Key 
tools to study this process through cosmic time are the galaxy 
stellar mass function (GSMF) and its integral over masses (the 
stellar mass density, SMD hereafter). 

Therefore, it is no surprise that most extragalactic surveys in 
the past decade have been used to determine the shape and evo- 
lution of the GSMF as a function of redshift. The earliest results 
based on small field surveys revealed that the SMD decreases 
with increasing redshift, as expected in the framework of the cur- 
rently accepted cosmolo gical hierarchical scenario, both in terms 
of integrated SMD (e.g | Giallongo et al.lll9 98: Dickinson et al. 
2003| iFontana et all 120031: iRudnick etalj 120031). as well as~of 
the de tailed GSMF (e.g. IFontana et al.ll2004t IDrorv et al.ll2004 
120051) . There have since been many indications that the evo- 
lution of the GSMF occurs more rapidly for more massive 



ysis to z > 3 (e. g. IFontana et alJ 120061: iKaiisawa et al.l [2009; 
Cap uti et al.ll2011l) . In parallel, wide-field sur veys have provided 
large samples with more accurate stat i stics dDrorv et al.l I200 " 1 



Pozz etti et all 1201 Ot IBolzonella et al.l l2010t iMarchesini et al. 
20 1 0: lllbertet al.ll2010ir One of the key results of these surveys 



galaxies than for low mass ones (e.g . Fontan aet al.1 200- 



Pozzetti et al.ll2007t Ipjrez-Gonzalez et al]l2008t IKaiisawa et al 
2009; IMarchesini et al.ll2.009t). a behaviour kno wn as downsiz- 



ing in stellar mass (see Fontanot et al.1 12009L and references 
therein). The advent of near- and mid-infrared facilities, above 
all the Spitzer telescope, has allowed the uncertainties in stellar 
mass estimates to be reduced and the extension of their anal- 
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has been the demonstration that the shape of the GSMF cannot 
be described by a (widely adopted) single Schechter function at 
least up to z — 1-5, but that it departs from this parametric form 
because of the superposition of individual distributions for the 
red an d blue galaxy pop ulatio ns dllbert et alJ2010l: | Pozzetti et alJ 
120101; IBolzonella et al.l l20"Toi; iMortlock et al.1 1201 lb or/and be- 
cause of a change with stellar ma ss either in star formation effi- 
ciency or galaxy assembly rate (IDrorv et alJl2009l) . 

An accurate knowledge of the GSMF is also a sensitive test 
of modern galaxy evolutionary models. From its initial stud- 
ies, much interest in the GSMF has been triggered by the pos- 
sibility of constraining the physics of the evolution of more 
massive galaxies, which, according to the hierarchical struc- 
ture formation scenario, are the results of the me rging of lower 
mass objects at earlier times dCole et alJ 1 1 9941) . In addition, 
to achieve a complete view of the galaxy formation picture, 
an important goal is a robust knowledge of the properties of 
low mass galaxies at high redshift. The slope of the GSMF at 
low masses may also represent a critical benchmark for current 
galaxy formation models. There is growing evidence that the 
number of low mass galaxies in the Universe is systematically 
overpredicted by most or all theoretical models dFontana et al.l 
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2006; iFontanot et all 120091: iMarchesini etai]l2009l) . Very sim- 
ilar eviden ce also appears in the analysis of the lum i nosity 
functions dPoli et alj 1200 1|: Marchesini & van Dokkuml l2007t 
lLo Faro et alJl2009HHenriques et al.ll201 lh . A particularly strik- 
ing aspect of the mismatch is that it appears in different ren- 
ditions of theoretical models, suggesting that it marks some 
fundamental incompleteness in our theoretical understanding of 
galaxy formation and growth. 

While a global picture is emerging from these investigations, 
many outstanding questions are still to be addressed. In gen- 
eral, the various GSMFs presented in the literature agree reason- 
ably well at z = 0-5, althou gh disagreements exist, somewhat 
increasing at high redshift ( Caputi et al. 201 it Gonzalez et alJ 



120 1 U IMarchesini et al. 2010; Mortloc k et alj|201 lh , that cannot 
be explained by merely field-to-field variance. At even higher 
redshift, the available estimates of the SMD are based on UV- 
selected samples, hence are potentially incomplete in mass, 
and/or are often derived by adopting average mass-to-light ratios 
for the w hole population rather than detailed object-by-object 
estimates dGonzalez et al.l [20 1 lh . Finally, and particularly rele- 
vant for the main topic of this paper, the GSMF at low masses is 
highly uncertain at intermediate and high redshifts, since current 
samples do not extend to the depths required to establish its slope 
with good accuracy. These uncertainties are due to a number of 
observational limitations. 

In addition to the uncertainties related to the GSMF compu- 
tation, it must not be forgotten that the actual estimates of stellar 
masses from broad band photometry are potentially affected by 
many systematic uncertainties, even when accurate redshifts are 
available. Part of this uncertainty is due to the lack of knowl- 
edge of important parameters of the stellar population, such as 
metallicity or extinction curve. The modelling of highly uncer- 
tain phases of stellar evolution is another source of uncertainty: 
in particular the different treatments of the thermally pulsating 
asymptotic giant branch (TP-AGB) phase is the source of the 
highe st discrepancies in simple stellar population models (see 
e.g. lMarastonll2005HMaripo et al.ll2008h . and has relevant impli- 
cations for the estimate of near-infrared luminosities and stel- 
lar masses for galaxies dominated by intermediate-age stellar 
populations (~ 1 Gyr). The largest bias is due to the difficul- 
ties in reconstructing the star formation history of each galaxy, 
which is necessary to estimate the appropriate M*/L ratio, and 
that may be poorly described by simplistic mode ls such as those 
adopted in stellar p opulation synthesis codes (lMaraston et alJ 
l2010HLee et aT1l2010h . 

All these uncertainties contribute to one of the main puz- 
zles that appear in present-day observational cosmology: the 
mismatch between the observed S MD and the integrated star 
formation rate density (SFRD) (e.g. | Hopkins & Beacoml [2006 ; 
iFardal et all 120071: IWilkins et al.ll2008l) . In principle, these two 
observables represent independent approaches to studying the 
mass assembly history from different points of view. However, 
the integrated star formation, after considering the gas recycle 
fraction into the interstellar medium, appears to be higher than 
the observed SMD at all redshifts. Several authors have high- 
lighted this severe d iscrepancy (of up to a factor of ~ 4 at z ~ 3, 
Wilkins et al. 2008). Moreover, if the mergin g contribution to 
the s tellar mass build-up is accounted for dDrorv & Alvarez 
2008), the agreement gets even worse. Intriguingly, the inte- 
grated SFRD exceeds the observed SMD, implying that we ei- 
ther overestimate the SFRD, or miss a substantial fraction of 
massive galaxies, or underestimate their masses, or finally fail in 
reconstructing the low-mass tail of the GSMF. An initial mass 
function (IMF) that varies over cosmic time was invoked to 



recon cile the two observables dFardal et al.ll2007t IWilkins et al.l 
2008). However, before invoking the non-universality of the 
IMF, it must be noted that both the SFRD and the SMD are 
affected by large uncertainties. The measure of the star forma- 
tion rate is itself particularly difficult, bein g either highly de- 
pendent on uncerta in dust corrections (e.g. ISantini et al J 12009: 
iNordon et al.ll2010h or limited to the brigh test far-infrared galax- 
ies at z < 2 - 3 (Rodig hiero et alj|20101) . On the SMD side, in 
addition to the uncertainties related to the stellar mass estimate 
itself, a major role is played by the poor knowledge of the low- 
mass tail of the GSMF. Owing to the limited depths of current 
IR surveys, the estimate of the faint-end slope basically relies 
on large extrapolations. An incorrect estimate, given the large 
number density of low mass objects, could translate into non- 
negligible errors in the SMD. 

A robust estimate of the slope of the GSMF is necessary to 
provide tighter constraints on all these unknowns. In this study 
we take advantage of the recent deep near-IR observations car- 
ried out by Wide Field Camera 3 (WFC3) installed on the HST in 
the upper part of the GOODS-S field in the Y, J and H bands and 
by Hawk-I mounted at VLT in the Kg band. These data allow ac- 
curate measurements of the stellar mass to very low limits. In this 
respect, we extend to hig her redshifts and lowe r masses the deep 

" d2009[) . The only study 
(2011), which was also 



analysis carried out by Kaiisawa et al 
of comparable depth is Mo rtlock et al. 



based on WFC3 data. However, the greater depth of the Early 
Release Science (ERS) images used in this work and the conser- 
vative cuts that we apply to the sample ensure an excellent over- 
all photometric quality, as we discuss in Sect. 13.31 Unfortunately, 
the area covered by ERS observations is small compared to re- 
cent surveys, and is slightly overdense. This feature somewhat 
limits the universal validity of our results regarding the SMD, 
especially in the intermediate redshift bins, although we chose 
our redshift intervals in order to ensure that the known clus- 
ters and groups (discussed in Sect. |2]i were mostly confined to 
two of them. However, we show that the study of the faint-end 
slope, which is the main aim of the present analysis, is insignif- 
icantly affected by these cosmic variance effects. In addition, 
this work represents an exercise to explore the potential of fu- 
ture deep WFC3 obs ervations, such as those of the CANDELS 
survey (Grogi n et al.ll201 it iKoekemoer et al. l l20Tll) . which will 
cover a much more extended area over various fields with depths 
comparable to the ERS observations. 

The paper is organized as follows: after introducing the data 
in Sect. [2] we present the stellar mass estimate and the GSMF in 
Sect. [3] the analysis of the faint-end slope in Sect. [4] the SMD 
and its comparison with the integrated SFRD in Sect. |5J and the 
comparison with theoretical predictions in Sect. [6j Sect. [7] sum- 
marizes our results. In the following, we adopt the A-CDM con- 
cordance cosmological model (Hq = 70 km/s/Mpc, Q.m = 0.3 
and Qa = 0.7). All magnitudes are in the AB system. 

2. The data sample 

This work exploits a new set of near-IR images that represent a 
significant improvement in photometric quality and depth over 
existing surveys. The first component is the public release of the 
Early Release Science (ERS) observations taken with WFC3, 
the new near-IR camera on board HST. The ERS observations 
cover an area of ~ 50 arcmin 2 , located in the northern ~ 30% of 
the GOODS-South field. They were taken in three filters, Fo98, 
7i25, and H^o, which reach 27.3 (Y) and 27.4 (J, H) magni- 
tudes at 5<x in an area of ~ 0. 1 1 arcsec 2 . We used the ERS mo- 
saics produced as described in IKoekemoer et ail d201 lh : we also 
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refer to Grazian et alJ (1201 lh for details of the catalogues, and 
IWindhorst et alJ (1201 ll) for a full description of the ERS obser- 
vational program. 

We complemented these images with new deep K$ band im- 
ages taken over the GOODS-S field with the near-IR VLT imager 
Hawk-I. The latter were taken in th e framework of a progra m de- 
signed to search for z ~ 7 galaxies (iCastellano et alJ2010allbl) . In 
the Ks band, the surveyed area covers 80% of the WFC3 ERS 
area. Owing to this and after excluding image edges of dubious 
quality, the available area reduces to ~ 33 arcmin 2 . The data re- 
duction of the Ks i mages is analogous to the procedure used for 
other Hawk-I data (Castella no et al.ll2010al) . The net exposure 
time is 25200s, with a lcr r.m.s. of 1.26 counts per second in a 
1" aperture. The magnitude limit at 5<x is ~ 25.5, one magnitude 
deeper than the previous ISAAC Kg band. 

We finally built a multiwavelength GOODS-ERS catalogue 
adding the other public images available in the GOODS- 
S field. They include the ACS images in the BVIz bands 
the deep UR images from VIMOS 
and the four IRAC bands at 3.6, 4.5, 



.. - 2004 
I2009h 



(Giavalisco et al. 
dNonino et al 

5.8, and 8.0 //m. With respect to the data set used to assem- 
ble our previous G OODS-MUSIC sample (IGrazian et al.ll2006t 
ISantini et al]|2009l) . the present GOODS-ERS data set benefits 
not only from the much deeper IR coverage provided by the new 
WFC3 and Ks band data, but also from a deeper version of the 
Z band image, which nearly doubles the exposure time of the 
previous image, a deeper U band image, and a brand new R im- 
age. In this data set, we extracted a 14 band multiwavelength 
catalogue using the H band as a detection image. Colours were 
carefully obtained with the same technique used in the GOODS- 
MUSIC catalo gue, where we adopt ed the PSF-matching code 
CONVPHOT foe Santis et alj|2007l) to accurately deblend ob- 
jects in the ground-based and Spitzer images. We note that the 
depth of the H band even exceeds the depth of the bluest bands, 
resulting in very poor quality photometric information about the 
faintest //-selected objects. 

The catalogue was cross-correlated with existing spectro- 
scopic samples. For sources lacking spectroscopic informa- 
tion, photometric redshifts were computed by fitting the 14 
band multiwavelength photometry to the PEGASE 2.0 templates 
( Fioc & Rocca-Volmerangd I1997L see details in IGrazian et alJ 
120061) . The accuracy reached by the photometric redshifts is very 
high, the absolute scatter |Az|/(l +z sp ec) being equal to 0.03, with 
only 3% of severe outliers (|Az|/(l + z sp ec) > 0.5). The statisti- 
cal error associated with each photometric redshift was used to 
evaluate the limiting magnitude at which a reliable GSMF can 
be computed. We found that a limit H 26 (or, equivalently, 
Kg — 25.5) is appropriate to maintain the error in the mass esti- 
mate (see next section) to within 0.3 dex and the relative scatter 
in the photometric redshifts |Az|/(l + z) < 0.1 for 85% of ob- 
jects, and we adopt this in the following. From the analysis of 
the individual photometric -redshift probability distributions, we 
can compute the fraction of "reliable" candidates. We considered 
a candidate to safely lie within a given redshift interval when the 
integral of its probability distribution curve, normalized to unity, 
over that interval is larger than 90%. Moreover, we accepted a 
certain level of tolerance in the definition of the redshift range 
to allow for the uncertainty in photometric redshifts. Following 
this method, for all Ks < 25.5 sources with z p hot > 2, we can 
exclude a secondary redshift solution at z p hot < 1.5 in 97.2% of 
the sources. This fraction increases to 99.6% when only bright 
sources (K s < 24) are considered. We also extracted a Ks band 
detection catalogue, and verified that all the objects detected in 
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Fig.l. Redshift distribution of the K s < 25.5 GOODS-ERS 
sample (black solid histogram) compared to t he Ks < 23.5 
GOODS-MUSIC one adopted bv lFontana elatl (12006!) (red dot- 
ted histogram). Overdensities at z — 0.7, z — 1, Z — 1.6, and 
z - 2.2 - 2.3 can be recognized in the distribution (see text for 
references). 



the Ks band are also detected in the H one, which is unsurprising 
given the extraordinary quality of the WFC3 data. 

On the basis of these results, we decided to restrict our analy- 
sis to the Ks < 25.5 sample, albeit obtained from the //-selected 
one, for two reasons: firstly, this selection allows a more robust 
comparison with previous ZT-selected surveys; and secondly, a 
Ks = 25.5 threshold is more efficient in detecting low mass 
objects than a H — 26 one. Adoptin g this cut, we exte nd by 
two magnitudes the previous work of iFontana et al.l (12006). who 
studied the GSMF of the K s < 23.5 GOODS-MUSIC sample. 
Our Ks < 25.5 sample here includes 3210 objects, 421 of which 
have spectroscopic redshifts. 

We plot in Fig. [T] the redshift distribution of the GOODS- 
ERS sample used in this work (black solid histogr am) compared 
to that of the GOODS-MUSIC sample adopted bv lFontana et al.l 
(2006) (red dotted histogram). Since the area covered by the ERS 
survey is relatively small, the sample is more sensitive to over- 
densities. The extended overdensities at z — 0-7 and z — 1, 
which cover the entire GOODS-S field (IVanzella et al.l [2005, 
ISalimbeni et"aT1l2009al and references therein), are clearly rec- 
ognizable. Unfortunately, the norther n part of GO ODS-S also 
includes a cluster at z - 1 .6 ([C astellano et al. 2007) a nd various 
groups at z 2.2 - 2.3 (ISalimbeni et al.ll2009at l Yang et alj2010t 
Magl iocchetti et al.ll2.01 ll) . which both affect the overall redshift 
distribution. 

Another difference from the Fontana et al. (2006) analysis is 
that the final sampl e used in this wor k includes Type 2 AGNs, 
since we show in Santin i et al ] (120121) that their stellar mass es- 
timate is insignificantly affected by the nuclear em ission. The 
same is untrue for Type 1 AGNs (Sant ini et alj |2012). so we re- 
moved spectroscopically identified Type 1 AGNs from the sam- 
ple. Since their number is very small (only four sources iden- 
tified in the entire sample), their removal does not affect the 
GSMF estimate. We also removed all identified Galactic stars. 
Finally, we applied a redshift selection to the range 0.6-4.5 and 
we ended up with a sample of 2709 objects (of which 354 have 
spectroscopic redshifts). 
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Fig. 2. Upper Distribution of the AM/M ratio, where AM is the 
average lcr error bar for each object, at z ~ 2 {left) and z ~ 3 
(right). Black solid histograms refer to the GOODS-ERS sam- 
ple, whereas red dotted ones represent the GOODS-MUSIC data 
set. Lower Median AM/M in each of the redshift bins used in 
this work as a function of the central redshift for GOODS-ERS 
(black solid circles/solid lines) and GOODS-MUSIC (red open 
boxes/dotted lines). 



3. The galaxy stellar mass function (GSMF) 

3.1. Stellar masses 

Stellar masses were estimated by fitting the 14 band photom- 
etry (up to the 5.5 fim rest-frame) to the Bruzual & Chariot 
syntheti c models in bo th their 2003 (BC03 hereafter) and 2007 
version (Bruzual 2007, CB07), through a^- 2 minimization. For 
consistency with our previous works and most of the studies 
in the literature, we adopted the estimates derived with BC03 
templates as the reference ones. In the fitting procedure, red- 
shifts were fixed to the spectroscopic or photometric ones. Our 
lcr errors, caused by both the photometric uncertainties and the 
photometric -redshift scatter, were computed by considering all 
the solutions within^- + 1 . During the error computation, spec- 
troscopic redshifts were fixed to their value, while photometric 
ones were allowed to vary around their best-fit solution in order 
to account for their degeneracy. 

We parametrized the star formation histories as exponen- 
tially declining laws with a timescale t. We analysed a wide pa- 
rameter space for metalli cities, ages, ex t inction s, and t, whose 
details can be foun d in Fontan a"et al.1 ([2004), as updated in 
Santi ni et al.1 d2009t) . With respect to our previous works, we 
also excluded templates with super-solar metallicity at z > 1 . 
Studies of the mass-metallicity relation dMaiolino et al.l 12008) 
indeed demonstrated that galaxies at high redshift are typically 
characterized by sub-solar metallicities. We decided to adopt ex- 
ponentially declining r models despite it being likely that they 
are a poor an d oversimplified desc ription for t he star formation 
history (e.g., iMaraston etall 12010b . However, iLee et all d2010l) 
showed that the resulting stellar masses can still be considered 



robust because of a combination of effects in the estimate of the 
galaxy star formation rates and ages. Moreover, r-models are 
widely used even in the most recent literature and allow a direct 
comparison with previous works. 

We adopted a Salpeter IMF. We also computed the stellar 
masses by assuming a Chabrier IMF, and checked that these 
are simply shifted by a factor -0.24 dex, which is constant to 
within 3% at the different redshifts. Moreover, we tested that the 
GSMFs obtained by adopting the two IMFs are consistent after 
applying the same shift , in agreement with what was found by 
Salim beni et"ai1 (l2009bl) . 

The comparison with our previous GOODS-MUSIC sam- 
ple allows us to test the effect of a deeper photometry data set 
on the accuracy of photometric redshifts and stellar masses. 
For this reason, we compared photometric redshifts and stel- 
lar masses for identical objects. The photometric redshifts of 
the present GOODS-ERS data are in very good agreement with 
the GOODS-MUSIC ones. Considering all objects in common 
between the two catalogues, the average scatter is < \zers ~ 
Zcoods-music\/(1 + Zers) >= 0.07, with only 0.06% of severe 
(scatter > 0.5) outliers. The stellar masses are also consistent 
with those derived from the GOODS-MUSIC catalogue. When 
selecting galaxies for which the redshift estimate differs by 0.1 
at most, the scatter (M ERS - M C00 ds-music)/Mers is on av- 
erage equal to -0.03 + 0.40. The major improvement provided 
by the higher quality photometry of WFC3 observations leads 
to a reduction in the uncertainties in the stellar masses. In Fig. 
[2] we compare the relative error in the stellar mass estimated 
using the GOODS-ERS data set (black solid circles/solid lines) 
with that obtained from GOODS-MUSIC catalogue (red open 
boxes/dotted lines). We again selected only galaxies common to 
both catalogues with consistent redshifts. In the upper panels, we 
show the distribution of the AM/M ratio, where AM is the aver- 
age lcr error bar for each object (AM = (M* „ mx - M* mi „)/2), in 
two redshift intervals centred at z = 2 and z = 3; in the lower 
panel, we have plotted, as a function of redshift, the median 
AM/M in each of the redshift bins used in this work. The relative 
errors in the stellar mass for the GOODS-ERS sample are on av- 
erage ~ 30% smaller than those computed for the same objects 
using the GOODS-MUSIC photometry, the difference increas- 
ing with redshift. It is clear that deep photometry in the near-IR 
regime is crucial to help improve our stellar mass estimates. 

3.2. The GSMF estimate 

We estimated the GSMF by ado pting both the non-parametric 
UVjnax method (ISchmidtl fl968l) and the STY dSandage et al.l 
1979) maximum-likelihood analysis assuming a Schechter para- 
metric form. As for any other magnitude-limited sample, our 
sample does not have a defined limit in stellar mass. For this 
reason, at each stellar mass and each redshift, we computed the 
fraction of objects lost because of the limited widt h of the M*/L 
distrib ution by adopting the technique described in lFontana et al.l 
(2004), after verifying that the simple parametrization used to 
describe the observed M*/L distribution still holds for our sam- 
ple. 

We show the results of our analysis as black solid circles 
(l/V max method) and black solid lines (STY approach) in Fig. 
[3] (where the reference BC03 templates were used). The best- 
fit Schechter parameters are reported in Tab. Q] The error bars 
in the \/V max points include Poissonian uncertainties, as well 
as uncertainties in the stellar masses. The latter were estimated 
by means of a Monte Carlo simulation, where we randomly ex- 
tracted the stellar masses according to their lcr uncertainties and 
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Fig. 3. GSMFs obtained with BC03 stellar templates in different redshift ranges compared with previous works. Black solid circles 

represent our analysis with the 1 /V max method, black solid lines show the best-fit to a Schechter function according to the STY 

approach. The grey dotted line replicates the best-fit Schechter function at z ~ 0.8 in the higher redshift panels. Error bars include 

the uncertainties in the stellar masses as well as Poissonian errors. The highest mass points are often poorly determined, because 

of the large statistical error, resulting from our poor sampling of the massive sid e. Other symbols rep resent the 1 /V max results 

of previous works, scaled to the same cosmology and converted to the same IMF: Fontana et al. (2006): open blue circles (F06); 

i — 1 j —j i 1 i —j i — ■ 1 i — i 

Bolzonella et al. (2010): solid green triangles (B10); Pozzetti et al. (2010): open purple triangles (P10); Ilbert et al. (2010): grey 

stars (TlO^ lKaiisawa et alJd2009l): orange cros ses (K091:lMarchesi ni et alTd2009l) : s olid red boxes (M09);LM archesi ni et alJ d2010): 

open pink boxes (M10); lMort lock et alT d201 ll) : solid dark green pentagons (Mol 1); Gonza lez et alJ d201 ll) : open violet pentagons 

(Gil). All the literature works considered for the comparison adopted the same stellar templates as this study. The legend shows 

the redshift intervals in which each set of points was computed. 



re-computed the GSMF of 10000 mock catalogues using the 
same procedure described above. 

Given the degeneracy between the faint-end slope a and the 
characteristic mass M* of the Schechter function, the STY ap- 
proach suffers from the incomplete sampling of the high mass 
regime owing to our small area, especially at high redshift. For 
this reason, the highest mass 1/V max points are often poorly de- 
termined, because of the large statistical error. In the highest red- 
shift bin, to constrain the fit, we M* to the value obtained 
at z ~ 3. The fits found in this case are shown in Fig.[3]and the 
relevant Schechter parameters are given in Tab. [T] In this table, 
we also provide (third column) the fraction of objects where a 



1 In Sect.f4] we study how the best-fit parameter a varies when choos- 
ing a different value of M* . 



secondary photometric -redshift solution falls outside each red- 
shift interval. This fraction was defined following the criterion 
discussed above and allowing a tolerance in photometric red- 
shift of 0.2. We found that the number density, given by the 
normalization parameter <p*, decreases with increasing redshift 
from 10" 3 - 70 -™ at z ~ 0.8 to lO -412 -™* at z ~ 4. We recall that 
cosmic variance effects could cause oscillations in the normal- 
ization parameter, especially in the two bins that are most af- 
fected by the presence of overdensities, namely the 1 .4 - 1 .8 and 
1 .8 -2.5 redshift intervals. However, a similarly decreasing trend 
for the normalization was also observed by prev i ous works (e.g. 
i Fontana et aLll2006t iPerez-Gonzalez et dj|2008t iKaiisawa et al.l 
120091: iMarchesini et al.ll2009HMortlock et al.ll201 ll) . Most iriter^ 
estingly, the low-mass slope steepens significantly from z ~ 0.8 
to z ~ 3, where the Schechter parameter a decreases from 
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-1.44+0.03 to -1.86+0.16, and then flattens fromz ~ 3 toz ~ 4. 
As demonstrated in Sect.|U this result remains valid here despite 
the uncertainties derived for the small area covered by our sur- 
vey and the presence of known overdensities. Indeed, even the 
redshift ranges that are the most contaminated show faint-end 
slopes in line with the results in the other redshift bins. 

Up to z ~ 2.5, the GSMF shows a dip at M* * 10 10 M o , 
which seems to shift to higher stellar masses as redshift in- 
creases, implying that a single Schechter is a poor parametriza- 
tion. This dip has been identified in previous wide-field sur- 
veys and interpreted as the differential evolution of the red 
and the blue populatio ns (Ilbert et al. 2010; Pozz etti et al.ll2010t 
iBolzonella et alj2010tlMortlock et al.l201 lb . The effect is larger 
in the redshift intervals 1.4 < z < 1.8 and 1.8 < z < 2.5, 
which are highly af fected by the presence of a well-known clus- 
ter at z ~ 1.6 (e.g. lCastellano et al.ll2007[) and of several of lo- 
calized overdensities at z — 2.2 - 2.3 (Salimbeni et aD|2009a; 
Yang et afll2010h . respectively: they are indeed populated by a 
higher fraction of old red galaxies, which enhances this dip. A 
di fferent explanation of the dip around ~ 10 10 M G was suggested 
by iDrorv et all (l2009h . who also measured a bimodal shape in 
the GSMF of the blue and red populations separately. This di- 
chotomy in galaxy formation, which pre-dates the red sequence 
appearance, was ascribed to a change with stellar mass in either 
star formation efficiency or galaxy assembly rate. The studies 
cited above show that a double Schechter is a more accurate de- 
scription of the shape of the total GSMF. However, given the 
small size of our sample, the inclusion of two more free param- 
eters (Bolzone lla et al.l l2010) makes the fit degenerate and was 
not an approach that we adopted here. 

3.3. Comparison with previous results 

We show in Fig. 0a compilation of 1/V max points collected from 
the literature, as listed in the legend, scaled to the same cos- 
mology and IMF. Unfortunately, it is impossible to correct for 
the effects of different stellar libraries, because, as we show in 
the next section, we are unable to determine a systematic shift 
that they could cause in the GSMF. Therefore, we decided to 
show only those studies that adopt the same stellar library as this 
work. Overall, our results are in good agreement with most of 
the other surveys, especially up to z ~ 3. In the two redshift 
intervals affected by the overdensities (z ~ 1.6 and z ~ 2.2), 
our GSMFs are on average higher than the other surveys, but 
still consistent with most of them within the errors. In general, 
we report a larger number of galaxies at the bright tail than 
iFontana et al.l d2006l) . because the present study includes AGNs 



in high mass galaxies (Bundv et al.l 2008 


; Alonso-Herrero et al.1 


2008; Brusa et al. 2009; Silverman et al. 


20091 Xue et all 20101 


Santini et al.l 2012b. Given the verv deep near-IR observations 



used in this work, the sampling of the low-mass end of the 
GSMF is considerably finer than most previous surveys, on av- 
erage by 0.5 dex up to z ~ 1.8 and by 0. 1 dex at z ~ 4, and at the 
same time the conservative photometric cut (K$ < 25.5) ensures 
reliable results even at the lowest masses. 

The only comparable study s ampling simila r or s lightly 
lower stellar masses is the one of Mortlock et aT| (1201 ll) . This 
work is somewhat peculiar, being obtained from a set of bi- 
ased pointings specifically designed to contain as many mas- 
sive galaxies as possible, and a posteriori corrected to account 
for this bias. They pushed their detection to H — 26.8 at a 5<x 
level, while our sample, although extracted from images of sim- 
ilar depth, was cut at a brighter limit to ensure good photometric 



quality. They also did not include any K band data, which is 
important to estimate reliable stellar masses. Finally, since our 
study is based on 14 ba nds of photometry (instead of 6 bands 
as Mortl ock et all 120 lib , our work also relies on good quality 
photometric redshifts. 

Despite the limited sky area, the bright-end tail is compara- 
ble overall within the uncertainties with that inferred by large 
surveys over the whole redshift range (with the exception of the 
1 .4 - 2.5 redshift interval, which, as discussed above, is affected 
by the presence of overdensities). The only severe disagreement 
is fou nd when comparing ou r results in the highest redshift inter- 
val to lGonzalez et al.l d201 lb . who, as already pointed out in the 
introduction, derived the GSMF by using a different procedure, 
i.e. by combining the UV luminosity function with an average 
M,/L ratio. 



3.4. The effect of different stellar templates 

The systematic uncertainties caused by the various assumptions 
involved in spectral energy distribution modelling were shown 
to dominate the overa ll error budget affecting the GSMF (see 
Marchesini et al. 2009, for a detailed analysis). In this regard, a 
significant role is played by the choice of the stellar templates 
used to estimate the stellar mass. 

Stellar masses obtained using the CB07 stellar library, which 
includes an improved TP-AGB stars treatment, are on average 
0.12 dex lower than those inferred using the BC03 templates, 
with a scatter as large as 0.17 dex. We plot in Fig. [4] their ratio 
as a function of the stellar mass adopted as a reference in this 
work (Msan) in different redshift bins. The lack of a clear trend 
of Msan/McBOi with stellar mass or redshift translates into a 
lack of a rigid offset between the GSMFs computed with the 
two libraries, although the CB07 points are on average at lower 
stellar masses than BC03. 

We compare in Fig. the GSMFs obtained with the BC03 
templates (black solid curves/solid circles) and the CB07 ones 
(red dotted curves/open boxes). For the sake of simplicity, we 
decided to report the four most representative bins. The results 
for the 1.0- 1 .4 and 1 .4 - 1.8 redshift bins are very similar to 



the 0.6 - 1.0 and 1. 8 - 2.5 ones, respective ly. We also show 
the l/V max points of iMarchesiniet all J2009h (their set 8) and 
ICaputiet aiJd2011b . b oth obtained by adopting the CB07 tern 



plates. The results of iMarchesini et all (2009) agree with our 
CB07 -based GSMF in all except the 1.8 - 2.5 redshift inter- 
val, likely because of im perfect redshift over lap between the two 
analysis. The points from Ca puti et"ai1 d201 lb are in broad agree- 
ment with ours at the bright end, while the incompleteness that 
the authors claim to be affected by below M, ~ 1O U M is likely 
responsible for the disagreement at low stellar masses. 

The best-fit Schechter parameters of the CB07 -based GSMF 
are reported in Tab. [2] At z > 2.5, we were forced to fi?<0 the M* 
parameter to its best-fit value at z ~ 2.15. If it is instead allowed 
to vary, the fit is unconstrained or the maximum-likelihood anal- 
ysis does not converge. The CB07- and BC03-based GSMFs dif- 
fer from each other. However, we do not find any similar sys- 
tematic behaviour at all redshifts. That the high-mass end of the 
CB07 -based GSMF is unconstrained at z > 2.5, while the BC03- 
based one suffers from poor statistical sampling only in the high- 
est redshift bin (z > 3.5), is a further confirmation that the two 
GSMFs are not affected by a systematic shift in stellar mass. 



2 See Sect. [4] for an analysis of how the best-fit parameter a varies 
when choosing a different value for M* . 
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Fig. 5 . Compariso n betw een the GSMFs obtained with the 
Bruzual & Charloll d2003l) template librar y (black solid circles 
and solid curves) and the iBruzuall (120071) one (red open boxes 
and dashed curves). Symbols are the results of the 1/V max anal- 
ysis and curves represent the STY Schechter fits. Error bars in- 
clude the uncertainties in the stellar masses as well as Poissonian 
errors. Other symbols present 1 /V max results of previous works 
based on CB07 templates, scaled to the same cosm ology and 
converted to the sa me IM F: Marc hesini et al.l (120091) : solid pur- 
ple triangles (M09): ICaputi et al.l d201 lb : cyan stars. 



At the lowest and the highest redshifts, we find the closer 
agreement, the normalization of the best-fit Schechter function 
being only slightly lower when CB07 templates are used. At in- 
termediate redshifts, we observe a more serious disagreement, 
resulting in different faint-end slopes and characteristic masses. 
This is unsurprising, because the effect of including of the TP- 
AGB phase is expected to be important at intermediate ages 
0.2 - 2 Gyr), which predominate the 2 < z < 3 redshift range 
Maraston 2005; Henriq ues et alj20"TH) . Although the difference 
between the CB07- and the BC03-based GSMFs do not show a 
systematic trend at all redshifts, the characteristic masses seem 
to be on average lower when CB07 templates are used (see Fig. 
[6] discussed in the next section), as expected, despite the large 
uncert ainties, in agreement with the results of iMarchesini et all 
(2009) . This trend is clear at 1.4 < z < 2.5, where o ur red- 
shift bins overlap with those of Marc hesini et al.l d2009l) . while 
the lack of high quality statistics at higher redshifts prevents us 
from drawing any firm conclusions about the effect of chang- 
ing the stellar templates. For what concerns the variation in a 
when changing the template library, we find similar slopes from 
z ~ 0.8 to z ~ 1.2, while in the redshift interval 1.4 - 2.5 the 
BC03-based GSM Fs are steeper t han the CB07-based ones by 
0.2 - 0.3. Marchesi ni et al.l ([2009) reports similar slopes when 
using BC03 and CB07 stellar templates in the redshift interval 
1.3 - 3.0, while their BC03-based GSMF is steeper than the 
CB07-based one at 3.0 < z < 4.0. The intrinsic difference be- 
tween the two surveys does not allow us to investigate the origin 
of this mismatch. 



4. The faint-end slope 

The main goal of this study has been to investigate the faint-end 
slope of the GSMF, especially at the highest redshifts (z > 2). 
From both Fig. [3] and Tables Q] and [2] it is evident that the low- 
mass tail steepens with redshift. The results from applying the 
STY approach to our BC03-based data indicate that the faint-end 
slope steepens significantly between z ~ 0.8, where we fitted 
a - -1.44 + 0.03, and z ~ 3, where the best-fit a is equal to 
-1.86 + 0.16, before flattening up to z ~ 4. 

First of all, we performed a simple sanity check to verify that 
the abundance of low mass objects at z > 1 .8 is reliable by plot- 
ting all objects with M* < 1O IO M and 1.8 < z < 2.5 on a BzK 
diagram. For galaxies at z > 2.5, we adopted the analogous RJL 
diagram (using IRAC 3.6 /-im as L band) , which extends th e for- 
mer to the 2.5 < z < 4 redshift regime dDaddi et al.l l2004). and 
checked stellar masses below 2 • 10 I0 M o . Approximately 91% 
of the sources indeed lie in the high redshift regions of these 
diagrams, making us confident of their photometric redshift esti- 
mate. As an additional check, we carefully inspected the individ- 
ual photometric -redshift probability distribution curves for each 
source with z > 1 .8 and M* < 2 ■ 10 10 M o . Following the criterion 
described in Sect. [2] we found that 96.5% of these sources have 
a 90% probability of lying at z > 1.5. 

As already pointed out in Sect. 13.21 the small sky area sam- 
pled by our data may be responsible for degeneracies between 
the faint-end slope a and the characteristic mass M* when fit- 
ting a Schechter function. We therefore studied in detail the de- 
generacies in the a - M* plane. The results are shown in Figs. [6] 
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Fig. 6. The a and M* space (lcr and 2cr contours) resulting 
from the maximum-likelihood analysis. Black solid curves refer 
to BC03-based GSMFs, red dotted curves refer to CB07 -based 
ones. 

andQ In the first figure, we analysed the redshift intervals where 
both parameters were allowed to vary, while in the second one 
we studied the dependence of the best-fit a on the chosen M* in 
those redshift bins where we were forced to fix the characteristic 
mass to constrain the maximum-likelihood analysis. 

In Fig. [6] we show the lcr and 2<x contours for a and 
M* Schechter parameters, for both the BC03-based GSMFs 
(black solid curves) and CB07-based ones (red dotted curves). 
While the parameter a is well-constrained at all redshifts (al- 
though with uncertainties increasing with z), our data prevent 
us from properly inferring the value of the characteristic mass. 
Nonetheless, as we show below, the result on a is robust against 
the degeneracy of M* . 

The steepening in a between z ~ 0.8 and z ~ 3 is clear 
from Fig. [6] when the BC03 stellar library is used. When we 
instead adopted CB07 templates, the faint-end slope did not 
change much from z ~ 0.8 to z ~ 2.2, while at higher redshifts 
(2.5 < z < 3.5) we were forced to fix the value of M* to constrain 
the fit (see Sect. I3.4l i. making the best-fit a parameter dependent 
on the choice of the characteristic mass. 

To study how the best-fit value of a changes when varying 
the M* value, we built a grid of log(M*[M ]) ranging from 10.5 
to 1 1.6 (these limits are justified by previous works) with steps 
in mass of 0. 1 and we fitted a Schechter function to the data for 
each element of the grid. We adopted this procedure in all the 
redshifts bins where the fit is unconstrained. We show in Fig. [7] 
the a - M* plane, where the shaded areas show the values of the 
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M* 

Fig. 7. Best-fit M* and a model parameters obtained by adopt- 
ing BC03 (upper panel) and CB07 (lower panel) templates at 
different redshifts. The symbol size increases with redshift. The 
shaded areas show the values of the faint-end slope at different 
fixed M* in the redshift bins indicated by the labels. 



faint-end slope at different fixed M* in the redshift bins indicated 
by the labels. The symbols show the best-fit values (and their un- 
certainties) for a and M* at z < 3.5 (z < 2.5 when using CB07 
templates), where our results are insignificantly affected by the 
lack of high quality statistics and we could allow both parame- 
ters to vary. The upper panel refers to the BC03-based GSMFs, 
while the lower one is obtained by adopting CB07 stellar tem- 
plates. 

From Fig(7](lower panel, blue shaded region), it appears that, 
whatever reasonable value for M* is chosen at 2.5 < z < 3.5, 
the best-fit a is clearly steeper than the best-fit values at lower 
redshifts, confirming the result found with BC03 templates and 
supporting that the major result of this paper is unaffected by 
the poor sampling of the high mass regime. However, as shown 
by the red shaded regions in Fig. [7] presently available data do 
not allow us to draw any firm conclusion about the trend of a 
between z ~ 3 and z ~ 4. 

To corroborate our result for the slope of the low-mass end, 
and ensure that is unaffected by the lack of statistics at the mas- 
sive end, we took advantage of the outcomes of large surveys 
and followed different approaches. We fitted the 1/V mnx points 
from this study together with those collected from the literature 
in comparable redshift intervals. We note that in principle, and 
in contrast to the STY approach, fitting l/V maA points involves 
data binning, thus may in general produce a different fit. We in- 
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Fig. 8. Faint-end slope as a function of redshift. The parameter 
a was computed through a maximum-likelihood analysis with 
a Schechter form (black solid circles refer to the BC03 library, 
red open boxes refer to the CB07 one) and by fitting the en- 
semble (this study + previous surveys) of 1/V max points with 
a Schechter parametric form (blue open triangles) and a double 
power-law shape (see text, green open circles). The different sets 
are shifted in redshift with respect to the central values in each 
interval (shown by the solid black circles). 



eluded only those surveys whose results are obtained using a 
method similar to our own and that sample the high-mass tail of 
the distribution, typically above M„ 3 x 10 10 M o . However, we 
obtained very similar results when also including the points from 
the literature at lower masses. We found that a single Schechter 
function does not seem to reproduce the faint- and the bright- 
end simultaneously in a satisfactory way. This is unsurprising 
because the Schechter function is itself a poor description of the 
shape of the GSMF when samples with high quality statistics 
are used (see discussion in the introduction and in Sect. I3.2l i. 
However, the inhomogeneity of the data set can also play a role: 
we collected 1 /V ma % points from different surveys, observed in 
different sky areas, and computed with slightly different meth- 
ods. We then fitted the ensemble of the l/V max points from 
this work plus those collected from the literature with a double 
power-lavvQ. The best-fit parameters are shown in Tab. [3] This 
analytic shape, having one additional degree of freedom than a 
single Schechter function, provides a tighter fit to the data at all 
redshifts. 

We report the different values of the faint-end slope as a func- 
tion of redshift in Fig. [8] It is shown that, regardless of the stellar 
templates and method adopted and the functional shape fitted to 
the data, all the results indicate a steepening of the faint-end of 
the GSMF with redshift up to z ~ 3. The trend is robust despite 
the relatively large error bars, especially at high redshift, and the 
presence of known overdensities at z ~ 1.6 and z ~ 2.2 - 2.3, 
and it is unaffected by the lack of high quality statistics at the 
high-mass end typical of small sky areas. The steepening of the 
faint-end slope with redshift seems to halt at z > 3 and the value 
of a seems to remain constant up to z ~ 4. However, although 
this is confirmed by the use of the outcome of previous large sur- 
veys, the results based on our data alone are largely dependent 
on the choice of the fixed M* parameter. 



3 The assumed functional shape is (/>,/[(M/M')- {,+a) + (M/MT <1+,8> ] 



The tendency for the low-mass end of the GSMF to steepen 
with redshift was previously fo und by other authors. According 
to the evolutionary STY fit of iFontana et al.l ([2006), a ranges 
from - 1 .25 ± 0.03 at z ~ 0.8 to - 1 .5 1 ± 0. 1 3 at z ~ 4. Increasing 
with redsh ift but flatter values for the faint-end slope were ob- 
tained by Mar chesini et al.l d2009f) , who reported a = -0-99 a t 
z ~ 1.6 and a = -1.39 at z ~ 3.5. iKaiisawa et all (2009) 
found a = -1.26+°;^ at z ~ 0.75 and a = -1.75+°;}f at 
z ~ 3. Very steep GSMFs, in agreement with those inferred 
with our data, were fitted by Capu ti et"aD d201 1[). who measu red 
a = -2.07+™ 8 , at z ~ 3.9 Finally. IMortlock et all d20Tlt) re- 
ported a = -1. 36 + 0.05 atz ~ 1.25 and a = -1.89 + 0.11 at 
z ~ 2.75. Both IMortlock et al.l (1201 ll) and lCaputi etafl (l201lh 
found a flattening at higher redshift similar to the one that we 
measure at z ~ 4, i.e. they fitted a - -1.73 ± 0.09 at z ~ 3.25 
and a = -1.85*g?2 at z ~ 4.6, respectively. 

We note that equally robust results cannot be inferred for the 
evolution of the characteristic mass M*, whose best-fit values are 
highly sensitive to the stellar templates and the functional shape 
fitted to the data (see Tabs. [T]|2] and [3}, as well as the size of the 
sample. 

5. The stellar mass density (SMD) 

We computed the total SMD by integrating the analytical fit- 
ting functions in each redshift bin from 10 8 to 1O 13 M . We 
show in Fig. [9] the SMD derived from the STY analysis using 
BC03 and CB07 templates as solid black and red curves, respec- 
tively. We also show a compilation of results from the litera- 
ture, reported to the same cosmology and IMF, as listed in the 
legend. The same integration limits as in this study were used 
in most of the works considered. The only excepti ons are the 
IMortlock etall d2011l) points (M, > 10 7 Mn). the lllbertetal. 
(20 IP]) one s (M„ > 10 5 M H ), and thos e from Ibickinso n et al. 
(2003|) and lPerez-Gonzalez et all (120081) . who adopted redshift- 
dependent mass limits (we refer to these works for more details). 
Our results show good agreement with those computed by pre- 
vious authors at 0.6 < z < 2, although we recall once again 
that our mass densities in the redshift intervals around z ~ 1 .6 
and z ~ 2.2 - 2.3 might be systematically too high owing to a 
few known overdensities. The steepness in the faint-end of the 
GSMF computed by this work is responsible for the large val- 
ues of the SMD inferred at z > 2. Our estimates are higher than 
those reported by pre vious authors, with the exception of the 
IMortlock et al.l d201 lb results. However, the latter results origi- 
nate from a different shape of the GSMF: IMortlock et al.l (120111) 
indeed found flatter faint-end slopes than we do, and the large 
SMD is a consequence of a higher density of high mass galaxies 
(see Fig.[3]l. 

We note that the SMD is affected by uncertainties caused by 
systematic effects. In Fig. [9] the grey shaded region indicates 
the dispersion in the SMD when including the outputs obtained 
by integrating the fit to our 1 /V max points plus those collected 
from the literature with both a Schechter function and a double 
power-law shape (see Sect. 0J. This region represents the sys- 
tematic errors caused by the choice of the stellar library and the 
functional shape of the GSMF, as well as the simultaneous use 
of the ERS observations as a probe of the low-mass end of the 
GSMF and the results of large surveys to constrain the bright- 
end. The dispersion increases significantly at z > 3, reflecting 
the large scatter among existing surveys. Moreover, the lack of 
overlap with most previous results at these redshifts is a sign 
of the impossibility to assemble a single, self-consistent GSMF 
from the highest to the lowest masses. This is due to the inhomo- 
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T QC +0.11 


8.23 


1.8- 


2.5 


425 


7.4% 


-1.84 ±0.06 


11.82 ± 0.28 


-4.17 +o-H 
-3.94 + ™ 


8.29 


2.5- 


3.5 


182 


11.8% 


-1.86 ±0.16 


11.30 ± 0.27 


7.97 


3.5- 


4.5 


51 


8.4% 


-1.80 ±0.20 


11.30 (fixed) 


a i T +8:38 

-0.10 


7.72 



Table 1. Best-fit Schechter parameters in the different redshift intervals as a result of the STY approach using IBruzual & CharlotI 
(120031) templates. Parameters with no error bars have been fixed to the value in the lower redshift bin. The second column indicates 
the numbers of galaxies in each redshift bin based on which the GSMF is actually computed. The third column shows the fraction 
of galaxies where a secondary redshift solution outside the redshift bin cannot be discarded with a 90% probability (see Sect. 13.21 
for details). The last column reports the corresponding mass density p obtained by integrating the GSMF between 10 8 and 10 13 M Q . 



Schechter parameters (STY method) - CB07 templates 



Redshift bin N %„„, secure a log 10 M*(M Q ) log 10 <MMpc 3 ) logp(M Q Mpc 3 ) 



0.6- 


1.0 


509 


8.5% 


-1.50 ±0.03 


11.85 ± 0.24 


-3.70 +01 ° 


8.39 


1.0- 


1.4 


372 


6.2% 


-1.46 ±0.05 


11.50 ± 0.18 


a 4 q +UI 

-Q 09 


8.22 


1.4- 


1.8 


264 


9.7% 


-1.42 ±0.07 


11.32 ± 0.17 


a 4| +0.00 
-Q \l 

Q CT +0.10 

"- ,,Ji -0.12 


8.09 


1.8- 


2.5 


437 


7.4% 


-1.58 ±0.06 


11.29 ± 0.15 


8.07 


2.5- 


3.5 


153 


11.8% 


-2.16 ±0.11 


11.29 (fixed) 


A qq +0.28 
^■ Jy -1.14 

-4 28 +ful 


8.05 


3.5- 


4.5 


45 


8.4% 


-1.88 ±0.21 


11.29 (fixed) 


7.67 



Table 2. Same as Tab. [T1using |Bruzua]| (120071) templates. 



Double power-law parameters (fit to 1/ V„ mx points) - BC03 templates 
Redshift bin a ft logi M*(M o ) log lo 0,»(Mpc~ 3 ) logp(M G Mpc~ 3 ) 



0.6-1.0 -1.36 ±0.02 -4.47 ±0.12 11.39 ±0.01 -2.75 8.50 

1.0-1.4 -1.52 ±0.02 -5.24 ±0.16 11.38 ±0.01 -3.10 8.24 

1.4- 1.8 -1.49 ±0.05 -4.70 ±0.23 11.30 ±0.03 - 3 - 12 -oofi 8 - 12 
1.8-2.5 -2.01 ± 0.04 -6.25 ± 1.57 11.64 ± 0.06 -3.94+™ 8.28 

2.5- 3.5 -2.28 ± 0.08 -6.70 ± 4.84 11.77 ±0.10 -4.73 8.24 
3.5-4.5 -2.27 ± 0.25 -6.38 ± 7.48 11.81 ±0.19 -4.84+™ 8.15 



Table 3. Best-fit parameters of the double power-law shape fit to 1 /V max points from this work (using BC03 templates) + a collection 
from the literature at M* > 3 x 10 I0 M o (see text). The last column reports the corresponding mass density p obtained by integrating 
the GSMF between 10 8 and 10 13 M . 



geneity of the samples, to the variance between different fields 
and also to the intrinsic uncertainties at high redshift in both the 
stellar masses and GSMF. 

We compared the SMD with the integrated star forma- 
tion rate density. For this purpose, we first considered the 
best-fit to the compilation of SFRD measurements mad e by 
iHopkins & Beacomld2006l) . Following Wilkins et alj d2008l) . we 
rescaled it to a Salpeter IMF and integrated it as a function 
of time, after accounting for the gas recycle fraction. The lat- 
ter is the fraction of stellar mass returned to the interstellar 
medium as a func t ion o f time, and was computed using the 
IBruzual & CharlotI (120031) model. The result of this calculation 
is shown in Fig. [9] by the blue dashed line. We then performed 
the same calculation by using the best -fit parametric shape fo r 
the star formation history inferred by Reddy & Steidel (2009), 
which also includes more recent high redshift points as well as 
a luminosity-dependent dust correction to the z > 2 data. We 
obtained the green long dashed line shown in Fig. [9] 

Our results solve the discrepancy between the SMD and the 
integrated SFRD at z > 2 (modulo the uncertainties affecting 



the z ~ 2.1 redshift interval), especially when considering the 
dispersion caused by the inclusion of high mass l/V max points 
fro m the other surveys. C onsi stency at high r e dshift was found 
by iMortlock et all d2011l) and iPapovich et all d2011l) . the latter 
study being based on an independents analysis. Overall, our re- 
sults support the notion that the SMD can be reasonably close to 
the integrated SFRD at z > 2, mostly due to a steepening of the 
GSMF, although our results might be systematically too high be- 
cause of the known overdensities in the small ERS field. As men- 
tioned, the higher values that we obtained than most previous 
studies is essentially due to the efficiency of WFC3 deep near- 
IR data to accurately recover the faint-end of the GSMF, espe- 
cially at high redshift, which contributes significantly to the total 
SMD. However, the significant steepening in the faint-end slope 
presented in this work is insufficient to solve the disagreement 
at z < 2, where the integrated SFRD exceeds the observed SMD 
by a factor of ~ 2-3, even when both of them are integrated 
down to low values of stellar mass / luminosity and adopting 
the SFRD computed by iReddv & Steidell d2009t) . Their SFRD, 
although it is lower than that resulting from the best-fit relation 
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Redshift 

Fig. 9. SMD between 10 8 and 1O 13 M versus redshift. The solid 
black and red curves are the integral of the STY Schechter fits 
using BC03 and CB07 templates, respectively. The grey shaded 
area shows the dispersion obtained when integrating the fit, with 
both a Schechter and a double power-law functional shape, of 
our l/y mflA points together with those from the literature (see 
text and Fig. [8]). Coloured symbols represent a compilation of 
results from the literature as listed in the legend (see text for 
their integration limits). All the points are scaled to the same 
cosmology and IMF. All the result s are b ased on BC03 stellar 
templates, except the Caput i et alj d2011|) w o rk, wh ich adopts 
the CB07 library. As far as iMarchesini et al.l d2010l) results are 
concerned, the central redshift was shifted by 0.02 for visualiza- 
tion purposes; given their Tab. 6, the central values is reported 
and the error bar indicates the total scatter in their estimates. 
The blue dashed and green long-dashed lines show the inte- 
grate d star formation history ac cordi ng to the best-fit r elatio n 
from lHopkins"& Beacom (2006) and iReddv & Steidell (f2009). 
respectively. 



of lHopk ins & Beacom (2006) at z < 2.5, is still unable to rec- 
oncile the two observables. The discrepancy is also not solved 
when our deep data, which allow a good control of the faint-end 
slope, are matched to the large surveys results to constrain the 
bright tail of the GSMF, and it gets even worse if one assumes 
that our mass density results are systematically too high owing 
to overdensities in the ERS field. 



6. Comparison with theoretical models 

In Fig. [T0l we compare our results with the predictions of semi- 
analytical models of galaxy formation and evolution, which fol- 
low the evolution of the baryonic component adopting an ap- 
proximate description of the relevant physical processes (i.e. 
gas cooling, star formation, stellar feedback, black hole growth, 
and AGN feedback) and of their interplay with gravitational 
processes, linked to the assembly of the large-scale structure 
of the Universe. These "recipes" include a number of param- 
eters that are usually fixed by comparing model predictions 



with a set of low-redshift observations. Despite their simpli- 
fied approach, semi-analytical models have turned into a flex- 
ible and widely used tool to explore a broad range of spe- 
cific physical assumptions, as well as the interplay between 
different physical processes. We considere d different, indepen- 
dently developed semi-analytical models: Menci etall 12 006) 
(red dotted curves) updated to i nclude thelReed et al.l (12007!) halo 
m ass function, MORG ANA dMonaco et al.l 120071, as updated 
Lo Faro et aD |2009, blue long -dashed curves), j Wang et al. 



(2008) (green dashed curves), and S omerville et alJ d2 01 ll) (or- 
ange dot-dashed curves). We refer to the original papers for a de- 
tailed description of the recipes adopted in t he galaxy formatio n 
and evolution modelling . All three models of lWang et al.l ([2008 ), 
ISomerville et alj d201ll) and MORGANA re solve galaxies with 
M« > 1O 9 M , while the iMenci et all (I2006T) model has a lower 
mass limit of 10 8 M o . The predicted stellar masses were con- 
volved with a Gaussian error distribution on logM, (see Fig. |2j 
to reproduce the observa tiona l uncertainties. W e refer the reader 
to iFontanot et al.1 d2009h and IMarchesini etal.1 (120091) for a de- 
tailed analysis of the effects of convolving the model predictions 
with observational errors. All stellar masses were converted to 
those for a Salpeter IMF. 

We first compared the GSMFs derived for various redshifts 
in the upper panel of Fig. |T0j we plotted only four redshift bins, 
the other two having very similar behaviours to the 1.0 - 1.4 
interval. We considered both the BC03- and the CB07 -based ob- 
served GSMFs. Despite the different physical recipes adopted 
by the different semi-analytical model s, their predictions ar e re- 
markably similar (as already found by Font anot et alJ2009h . All 
the models considered consistently predict a larger abundance 
of low mass galaxies than observations at least up to z ~ 3, de- 
spite the steep low-mass end slop e inferred from our data. The 
only exception is for the model of I Wang et al.l (|2008), which is 
closer than all other models to the observations at z ~ 2.2 and 
consistent with them at z ~ 3. The general overestimation of 
the faint-end is often attributed to a too efficient formation of 
low-to-interme diate mass (10 9 — IO'^q) galaxie s in the models 
dFontanot et al.ll2009l see also lLo Faro et al.ll2009l) and it cannot 
be explained by systematic uncertainties caused by the stellar 
templates. On the observational side, it is unlikely that so many 
low mass galaxies have been missed by observational surveys, 
especially at low redshift (z ~ 1), when the disagreement is ev- 
ident for M, > 10 1 M o . In our highe st redshift bin (3.5 < z < 
4.5), the models of IWang et all d2008l) . ISomerville et all d201 ll) 
and MOR GANA show a re asonable agreement with the data, 
while the Menci et al. (2006) model still slightly overpredicts the 
space density of low mass galaxies. This is because the model of 
IMenci et"al] J2006) includes starburst events triggered by fly-by 
interactions, which are very common in low mass objects at high 
redshift and increase their stellar mass. All theoretical predic- 
tions also underestimate the stellar mass of high mass galaxies 
in the highest redshift bin. However, this mass range is highly 
affected by cosmic variance, and the small (< 3) number of 
galaxies observed with our data in each of the highest mass bins 
(M > 10 n M p) prevents us from draw ing any firm conclusions. 
We note that Marchesin i et a 11 (120091) found a similar disagree- 
ment at high masses, using observations on a sky area ~ 10 times 
larger than the area sampled by this study. 

A different treatment for the TP-AGB phase in the stellar 
population synthesis models (i.e. moving from BC03 to CB07) 
may indeed help in reducing the high-mass end discrepancy, by 
redu cing the inferred stella r masses of the galaxies in the sam- 
ple. [Henriques etal] d201 ll) indeed demonstrated that the inclu- 
sion of a more detailed treatment of the TP-AGB phase in the 
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Redshift 



Fig. 10. Observed GSMFs in different redshift ranges (upper) 
and observed SMD between 10 9 and 1O 13 M (lower) compared 
with theoretic al predictions , show n as coloured curves. Red 
dotted curves: jMenci et alj d2006l). Blue long-dashed curves: 
MORGANA dMonaco et all l20Q7h . Green dashed curves: 
I Wang et all (12008). Orange dot-dashed curves: [Somerville et al. 
(120111) . Black solid circles, solid lines, red open boxes, and the 
grey shaded area show the results of the present work as in Figs. 
@and[5] 



semi-analytical framework solves a comparable mismatch be- 
tween the predicted and observed rest-frame K band luminosity 
function of bright galaxies at z ~ 2 - 3, by increasing the pre- 
dicted /^-band flux in model galaxies. This effect is particularly 
relevant at these cosmic ages, since model galaxies are domi- 
nated by intermediate-age stellar populations, with the right age 
(~ lGyr) to develop a TP-AGB phase, which makes them redder 
without the need to change their mass or age. We test this idea by 
considering stellar masses estimated using the CB07 synthetic 
templates. However, as clearly shown in Fig. [TUl (upper panel, 
red open boxes), assuming stellar templates that account for TP- 
AGB stars alleviates (the number of M > lO n M galaxies in 
the sample is reduced by 37%) but does not completely solve 
the mismatch at the high-mass-end of the GSMF in the highest 
redshift bin (where the previous percentage reduces to 12.5%). 

We also compared the various SMDs for galaxies of mass 
10 9 < M„/Mq < 10 13 , as predicted by the four models we con- 



sider (lower panel of Fig. fTTTb. and we found that it is on aver- 
age higher tha n that observed by a factor of ~ 2 up to z ~ 2. 
The models of IWang et all (120081) . ISomerville et al.1 (1201 ll) . and 
MORGANA show good agreement with our results at 2 < z < 3, 
and are below the data at the highest redshifts. The predic- 
tions of iMenci et al.l (120061) . given its overall steeper low-mass 
end, are instead higher than the observed SMD up to z ~ 3. 
These results are driven by the afore mentioned overabundance 
of intermediate-to-low-mass galaxies, which dominate the SMD 
at all redshifts, and are counterbalanced by the underestimation 
of the stellar mass in massive galaxies at high redshift. 

7. Summary 

We have used deep near-IR observations (Y as faint as 27.3, both 
J and H as faint as 27.4 AB mag at 5<x) carried out with the 
Wide Field Camera 3 in the GOODS-S field as part of the Early 
Release Science. These data, complemented with a deep K$ (as 
faint as 25.5 at 5<x) Hawk-I band data set and high quality pho- 
tometry in various bands from the near-UV to 8 jum, have been 
used to derive accurate estimates of the stellar mass. We have 
succeeded in reducing the average relative error in the stellar 
mass for a given object by ~ 30% with respect to that obtained 
with the GOODS-MUSIC catalogue in previous works from our 
group. 

Unfortunately, the sky region covered by ERS observations 
is biased b y a number of localized as well as diffuse overdensi- 
ties dVanzella et al.l2005|:ICastellano et alj2q07HSalimbeni et al.l 
l2009aHYang et al.ll2010l:lMaghocchetti et all201 ll) . which could 
be responsible for the oscillations in the normalization of the 
GSMF given the limited size of the area. However, this data set 
offers a unique combination of accuracy in photometric qual- 
ity (hence in photometric redshifts and stellar mass) and depth, 
which makes it ideally suited to studying the faint-end of the 
GSMF. 

We computed the GSMFs in six different redshift intervals 
between 0.6 and 4.5. Thanks to the depth of the catalogue, we 
were able to study the low-mass end of the GSMF at lower 
masses than most previous studies by 0.5 dex up to z ~ 1.8 and 
0. 1 dex at z ~ 4. We compared our results with previous works 
and found general good agreement, even at the highest masses, 
despite the limited sky area sampled by our data set. We found 
that the only redshift intervals that show poorer agreement with 
previous results are those between z = 1.4 and z = 2.5, which 
are known to be affected by the presence of overdensities: we 
discovered, as expected, a slightly larger abundance of massive 
galaxies. We also compared our results for the GSMF obtained 
with two different stellar libraries, BC03 templates and CB07 
ones, the latter including a treatment of TP-AGB stars. The stel- 
lar masses inferred from the CB07 library are on average 0.12 
dex lower than the BC03-based ones, with a large (0.17 dex) 
scatter. The lack of a clear trend with stellar mass or redshift of 
the ratio of the two estimates translates into a lack of systematic 
difference between the best-fit Schechter parameters in the two 
cases. The largest disagreement was found at z ~ 2-3, where the 
effect of the TP-AGB phase is expected to be the most important. 

The main result of this study is the steepening of the faint- 
end slope: the value of a increases from -1.44 ± 0.03 at z ~ 0.8 
to -1.86 + 0.16 at z ~ 3, and then flattens up to z ~ 4. We have 
confirmed the steepening of the low-mass end, which had been 
pointed out by previous authors, with deeper and higher quality 
photometry. Our results are unaffected by degeneracies in the 
M* parameter, and they are insensitive to the choice of either 
the stellar templates or the functional shape fitted to the GSMF, 
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as well as to the limitations of the small area covered by ERS 
observations. 

We computed the SMD as a function of redshift and com- 
pared it with the int e grated star formation histories derived by 
iHopkins & Beacoml (120061) and iReddv & Steidel d2009h . The 
finer sampling of the GSMF at low masses and the steep in- 
ferred faint-end slopes determine the higher SMD estimates at 
Z > 2 than most previous works, solving the disagreement ob- 
served by previous authors between the SMD and the integrated 
SFRD at these redshifts. However, despite the steep GSMF that 
we find, the integrated star formation history still exceeds the di- 
rect measure of the SMD at z ~ 2 by a factor of ~ 2 — 3, even 
when our data are analysed together with the results of previous 
large surveys to ensure a good sampling of also the bright-end 
tail of the GSMF. 

Finally, we compared our GSMF and SMD estimates with 
the predictions of four models of galaxy formation and evolu- 
tion. All models predict a larger abundance of low mass galaxies 
than observations, at least up to z ~ 3. They also underestimate 
the stellar mass of high mass galaxies in the highest redshift 
bin, although cosmic variance effects prevent us from drawing 
firm conclusions at these redshifts. The overabundance of low 
mass galaxies translates into a general overestimation of the to- 
tal SMD with respect to the data up to z ~ 2, while this density is 
underestimated at z > 3 owing to the dearth of massive galaxies 
at these redshifts. The exact degree of disagreement depends on 
the particular model. 

Future CANDELS data will cover a larger sky area and al- 
low a finer sampling of both the bright-end of the GSMF and its 
normalization, and at the same time they will be deep enough to 
accurately probe the GSMF faint-end. These, together with spec- 
troscopic follow-up campaigns, will reduce the uncertainties in 
the stellar masses, and they will significantly improve our results 
and our understanding of the stellar mass assembly process. 
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